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Noiseless subsystems offer a general and efficient method for protecting quantum information in the presence 
of noise that has symmetry properties. A paradigmatic class of error models displaying non-trivial symmetries 
emerges under collective noise behavior, which implies a permutationally-invariant interaction between the 
system and the environment. We describe experiments demonstrating the preservation of a bit of quantum 
information encoded in a three qubit noiseless subsystem for general collective noise. A complete set of input 
^Nj ■ states is used to determine the super-operator for the implemented one-qubit process and to confirm that the 

' fidelity of entanglement is improved for a large, non-commutative set of engineered errors. To date, this is the 

' largest set of error operators that has been successfully corrected for by any quantum code. 
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I. INTRODUCTION 

> 

f — Quantum information processing (QIP) holds the promise of solving problems in quantum simulation, quantum computation, 

I/"} and secure communication that have no known efficient solution in classical information processing While QIP can be, 
in principle, abstractly characterized without reference to the details of a specific implementation, physically realizing quantum 
information and its manipulation is essential to practically exploit its unique capabilities [^]. Real physical systems are invariably 
exposed to environmental noise and decoherence due to incomplete isolation from their surroundings, as well as to operational 
errors caused by imperfect manipulations. Thus, physical realizations of QIP are confronted with the challenge of achieving 
noise control during storage and processing of quantum information. 

Thanks to a series of recent investigations, major progress has been witnessed on the theory front of reliable QIP. On one 
side, powerful "accuracy threshold theorems" for fault-tolerant quantum error correction (QEC) [|[ 0. § 0] ensure that, if the 
^ integrated effect of the noise per qubit and computational step remains sufficiently small, arbitrarily accurate QIP is still possible, 
^ in principle, with reasonable resource overheads. On the other side, alternative noise control techniques have become available as 
^ options complementing and expanding the applicability of conventional QEC. These methods include both passive error control 
codes based on decoherence-free subspaces (DFSs) 0JlO| and noiseless subsystems (NSs) | pT| |l2| |l3|, |T^, and active error 
suppression schemes relying on dynamical decoupling [|15|,|16| , [T^ along with a variety of device-dependent schemes 

■ f-n 

X. 

H , is the realization that protecting quantum information against noise need not require the overall state of the physical device 
- - ■ supporting QIP to be perfect [|5|. Defining in what sense quantum information can be accurately stored in the noisy state of a 
physical system leads to think of all possible noise control options in terms of the emergence of logical subsystems (or "abstract" 
quantum particles) which are or can be made immune to noise |l2| ^ . In the former case, the occurrence of a NS directly 
ties into the existence of symmetry properties of the natural noise process, and error-free information storage is ensured without 
requiring active intervention |^ |ll]]. In the latter case, an external control action is instead necessary to appropriately symmetrize 
the dynamics [|l2, 14, 2^] or to enforce noiselessness by a process that incorporates active recovery operations as well [|ll]|. The 



for reducing systematic and calibration errors ^ ^ 

In retrospect, what constitutes the unifying conceptual feature of these advances, and what ultimately made them possible. 



basic intuition remains, nevertheless, unchanged. 

While the implications of the subsystem approach are still being investigated, the idea of separating the abstract information- 
carrying degrees of freedom from the implemented qubit degrees of freedom has proven useful beyond the original motivation 
of gaining noise protection. In particular, the notion of a subsystem has led to an operational prescription for realizing logical 
qubits in physical systems |^], and the analysis of the resources required to universally control subsystem-encoded qubits has 
resulted in an encoded universality approach [|2| ^8|], where the notion of universality may be tailored to the set of physically 
available interactions. More recently, the idea of a NS has shown to also underlie topological approaches to QIP | 
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From the experimental point of view, the significance of the NS notion has only recently begun to be explored. In particular, 
the first implementation of a non-trivial NS for general collective noise in a three-qubit liquid-state nuclear magnetic resonance 
(NMR) quantum information processor was reported by Viola et al. [Q. Here, we deepen our earlier investigation in two ways: 
by presenting an expanded description of the theoretical analysis and experimental methodology underlying [^o|]; by reporting 
additional experimental results which may further shed light on the relevance of NSs within present-day quantum information 
technology. 



II. THEORETICAL BACKGROUND 



A. Noiseless subsystems: From simple examples to the general definition 



We begin by trying to build the intuition underlying the NS idea based on simple considerations and prototype examples 
1^, It is a well-known lesson in physics that the occurrence of symmetries in a system generally implies the existence of 
"conserved quantities", and that these can be exploited to ease the understanding of the system's behavior In its essence, the NS 
approach adapts this lesson to the QlP-motivated task of achieving protection against noise. This is done by noticing that even 
though the system as a whole may be experiencing errors, some of its properties may still remain unaffected by them. Thus, if 
information can be represented in terms of the abstract degrees of freedom (DOFs) corresponding to such conserved quantities, 
noiselessness may be ensured in spite of the fact that the errors evolve the overall system's state. 

This intuition applies to both classical and quantum information storage. A simple classical example arises by considering 
two bits subject to errors which either flip both bits with probability p, or leave them alone with probability I — p. This error 
model has the property that it preserves the parity P{s) of a state s of the bits, where P{s) is defined as the sum (mod 2) of the 
bit string s. Thus, the two values of P can be used together with the (non-conserved) value of, say, the first physical bit to label 
the four possible states of the two bits: 

00 ^ 0-0 

01 ^ 1-0 

10 ^ 1 • 1 

11 ^ • 1 . 



The above table establishes a correspondence between the state space of the physical system (left hand side) and the one of a 
pair of abstract subsystems (right hand side). In the resulting representation, the first member of the abstract pair (the parity bit) 
carries the information to be protected, while the second (the "syndrome" bit) experiences the effect of the errors. In this case, 
parity provides a classical NS, and the information resides in this protected DOF. 

Many of the above features carry over to the quantum-mechanical case. In quantum systems, the presence of symmetries 
(hence of conserved quantities) is associated with the existence of operators that commute with all possible errors. First, consider 
a simple two-qubit example which corresponds to complete depolarization on qubit 2. This error model is defined by the set of 
error operators {Ea} = {-5 cr^\ a = 0,x, y, z}, aj^^ = Clearly, no joint state of the two qubits is protected. However, a 
one-qubit state can be stored in the first physical qubit without being affected by the errors: qubit 1 is a "trivial" quantum NS. 
Mathematically, this intuition is made precise by observing that, because [Ea, cri^''] = for every error and u = x,y, z, all the 
expectations of tri^^ - hence the state of qubit 1 - are protected from noise: ai^' ,u — x,y, z, define the observables of this trivial 
NS qubit. 

In general, knowing the symmetries of the error model suffices for identifying possible protected DOFs, in a way similar 
to what happens in the above simple example: one first determines the set of operators that commute with the possible errors, 
and then in this set identifies appropriate combinations that algebraically behave like the observables (the Pauli operators) for 
abstract qubits If S" is a (finite-dimensional) quantum system interacting with some environment E, a description of noise on 
S which suffices for discussing error protection and error correction properties can be obtained by constructing the appropriate 
interaction algebra A JTl]]. Let the open-system evolution of S involve coupling operators Ja, where the Ja are traceless and we 
assume, for the moment, that the internal Hamiltonian Hs of S alone can be either set to zero or identified with one of the Ja- 
Then A contains all the complex linear combinations of arbitrary products of the Ja and the identity. If S is initialized in the 
state gin, and the evolution is depicted in terms of a quantum operation 

Qin '-^ Bout = £{Bin) = ^ EaQinEl , ^ E^Ea = 1 , (1) 
a a 

then A collects all the possible errors that the coupling to E can induce for arbitrary strength or interaction time. The commutant 
A' of A, which collects all operators on S that commute with arbitrary errors in A, is the relevant structure to be used in 
generalizing the symmetry argument given above. By construction, both A and A' are (multiplicative) sub-algebras of the full 
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algebra As of operators on S. Within the Hamiltonian framework for the composite system S, E adopted here, they are also 
naturally closed under the f operation, making some standard results from the representation theory of operator algebras directly 
applicable ^ Q|. The symmetry properties of error models corresponding to A can be visualized by thinking of the largest 
group of unitary operators contained in A' as the symmetry group for the problem. 

Suppose that S consists of n physical qubits, in which case the state space Ti. ~ (C^)®", and As can be identified with the 
algebra of complex matrices Mat7v(C) acting on C^, N = 2'"-. Then three possibilities are worth examining: 



If A' consists only of scalar multiples of the 1, no useful symmetries are present. Mathematically, this condition is 
equivalent to the fact that the error algebra A acts "irreducibly" on H, hence A ~ As by the Schur lemma JTl] , |3^ . 
Because the error process is contributed by all operators on S*, no protected DOF exists. 

Suppose instead that A! is non-trivial, which implies that ^ is a proper sub-algebra of As = Matjv(C). If ^ ~ Mat„i(C) 
for some integer to, then one can show [ ^ ] that there exists a change of basis U on 7Y, and integers d, r with md + r ~ N, 

U: C" ® C'' ® C^ (2) 

such that, in the resulting representation, all error operators leave the factor C™ unaffected: 

UAU^ = hn ® Matd(C) e Or , (3) 

where Or accounts for the action on the remaining summand C". Situations where a NS can be directly identified with 
a subset of the physical qubits belong to this category. The above trivial two-qubit example, for instance, corresponds to 

TV = 4, m = d = 2,r = 0. 

If A is again strictly contained in As, but it cannot be identified with the operator algebra of a m-dimensional quantum 
DOF, a theorem from the representation theory of operator algebras (the Wedderburn theorem, [ p3| , ^4| ]; see also [^) 
still implies that a transformation [/ to a new basis exists where, as above, the action of the errors takes a simple block- 
factorized form: 

[/ : -> C™i ® . . . ® C"^ , (4) 

with J2j ''^jd-j — and 

UAU'' ^Irm® Matdi (C) ® . . . ® 1™, ® Matd, . (5) 

Because, in this representation, the action of any error operator in A is only experienced by the "syndrome factors" 

in the decomposition (Q), each of the co-factors C™^ can be identified with the state space of a to^ -dimensional NS under 



It is worth noting that, with respect to the same basis (Q), the action of operators commuting with the errors becomes "dual" 
to the one of the errors themselves, 

UA'U'' = Mat™, (C) ® ® . . . ® Mat™, (C) ®td, , (6) 

implying that a non-trivial transformation is now enforced on the noiseless DOFs. Also, whenever dj = 1 for some j, the 
syndrome subsystem becomes effectively a classical DOF with a one-point configuration space. Thus, C™^ (K) C ~ C™^ , and 
the jth summand in the decomposition (Q) can accordingly be identified with a DFS under A [|ri], |l4|. 

Because of the abstract algebraic nature of the subsystem identification given by the mapping between the states of 

the information-carrying NSs and the states of the underlying physical qubits may in general become very indirect. However, 



the method for constructing a NS from appropriate observables in A' may be applied in general 126|. We now make these 
considerations explicit in the situation that is relevant to the experimental implementation. 



B. Collective noise for three qubits 

The collective error behavior provides the paradigmatic situation for discussing passive noise control through both DFSs and 
NSs, and in particular for realizing the simplest non-trivial noiseless quantum subsystem. Collective error models have been 
extensively analyzed in the theory literature [HJlO JIT 2^, and experimentally investigated in the QIP context in optical 
[p^], trapped-ion [^, and liquid-state NMR [pO|, [40|] devices. For a system S composed of n qubits as above, collective noise 
behavior arises whenever a single environment E couples to the individual particles without distinguishing among them. This 
results in error models that are characterized by permutation symmetry. Whether or not the natural dynamics of S alone, ruled 
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by Hs, actually respects this symmetry (as assumed so far) is an important issue with both conceptual and practical implications. 
While deferring a more detailed discussion of this point to a later stage, we begin by examining the consequences of permutation 
symmetry. 

In particular, let us focus on the situation where S is composed of three qubits, implying that H ~ C^, and the full operator 
algebra ^5 ~ Mat8(C). For the purpose of characterizing collective error models, the crucial property is that only global error 

generators J„ = (crl^'' + cri^'' + a';^^)/2, u — x, y, z, may be present in the system-environment interaction. By definition, Ju 
is the projection of the total spin angular momentum along the u-axis (in units h). In NMR with spin- 1/2 nuclei, for instance, 
interactions of this type may arise from uniform, fully-correlated magnetic fields which fluctuate in direction and strength, 
leading to evolutions which can be semi-classically described as 

I*)i23 ^ |*'>i23 = e-^'^'-'^+'y'y+'-^'^\^),2^ , (7) 

for random variables Ou, u = x,y, z. In the formalism of quantum operations, collective error processes are characterized by 
completely positive dynamical maps of the form (jl]), where the possible errors Ea are constrained to commute with all possible 
particle permutations. Thus, the largest interaction algebra Ac resulting from arbitrary collective interactions consists of all the 
totally symmetric operators on three qubits. Because the dimension of the subspace of totally symmetric operators for n qubits 
is given by {n + l){n + 2){n + 3)/6 |ft2|], Ac is a 20-dimensional sub-algebra of ^^5. This implies that a description of the most 
general collective error model on three qubits can be accomplished by using an error basis with at most 20 (linearly independent) 
operators, out of the possible 64 needed for representing arbitrary noise in the absence of symmetries. By writing 

Ac = span{E„ | a = 0, . . . , 19} C y^s , (8) 

an explicit basis of operators Ea can be constructed by fully symmetrizing the standard Pauli product operator basis for three 
qubits. Let us introduce compact notations to describe operators that are invariant under the full set of qubit permutations: 

. a^^a!^^ + .(^)a(^) + + ,(i),(3) + ,(2),(3) + ,(.),(3) , 

and so forth. Then a basis for Ac is given by the 1, the three linear operators J^, the six quadratic operators XX, YY, 
ZZ, XZ, XY , Y Z, and the ten cubic operators resulting from the above construction. 

As noted earlier, the symmetry properties of a given error model appear explicitly in the commutant of the error algebra. In 
the case of Ac, the commutant A'c contains the subgroup 11 of unitary operators that implement permutations of the particles 
e.g., a swap 7ri2 between qubit 1, 2 means 7ri2|«)i|j)2|fc)3 = \i) 2\j) i\k) 71, — |j)i|i)2|fc)3- In fact, one can show that the whole 
A'^ consists of linear combinations of operators in 11, expressing the fact that A'^ coincides with the group algebra CiSa of the 
permutation group ^3 under the above representation in H. Jl4| , |2^ , |3^ ] . Physically, because the error generators Ju are also the 
generators for the global rotations of the qubits, A'c can be regarded as containing the operators which remain invariant under 
such rotations. Thus, operators in A'c can be constructed from the identity and the simplest invariant operators [|[ pi] ]: 

S12 = ^(1) -at^), S23 = O^(') • S31 = ^(3) • <7(l^ (10) 

where • denotes the usual dot product. Note that Sjk is nothing but the Heisenberg spin coupling between spins j, k. As it turns 
out, every operator in A'c can be realized, in principle, through the application of Heisenberg Hamiltonians of the form (|l0|) [|43|]. 
By observing that the total angular momentum observable ,P = J ■ J simply rewrites in terms of the above operators, collective 
symmetry immediately implies a conserved quantity, given by the eigenvalues j (j + 1) of J^. 




C. Abelian error models 

Suppose that the system-environment interaction is contributed by a single global error generator = J ■ v, v ■ v — 1. Then, 
in the absence of an independent quantizing direction (provided, for instance by Hg), the resulting interaction algebra is abelian, 
and the corresponding error models accounts for collective decoherence (unrecoverable loss of phase information) with respect 
to the fixed basis of eigenstates of (jy. The choice v — z singles out the z basis. For three qubits, the relevant error algebra Az 
can be constructed from the generator and the 1, and can be identified with the sub-algebra of Ac spanned by the four axially 
symmetric and permutation-invariant operators, i.e. 



Az = 5pw{t,Jz,'ZZ,ZZZ} c Ac , 



(11) 
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For instance, an error model in this class that will be of practical significance is a full-strength (or "crusher", borrowing from 
the NMR terminology) collective z-dephasing on three qubits, which may be described by a quantum operation with Kraus 
operation elements K^, a = 0, . . . , 3: 



= i (l + 2 + + ZZZ^ , 
Kl = i (l - 2 + - ZZZ] , 



X| i + 2J^-'ZZ - iZZZ^ , 

= ^{?,l~2J^-ZZ + 2,ZZZ^ . (12) 

In a similar way, collective decoherence about the x axis can be described by effectively switching from the z to the x basis, 
i.e. by mapping Az ^ Ax = HAzH, where H represents a collective Hadamard transform. Accordingly, a crusher collective 
x-error process corresponds to = {-^f}- with operation elements obtained from ( p^ via the appropriate rotations. 

For later purposes of comparison between encoded and un-encoded information, the description of the noise process induced 
by a three-qubit error model on a physical information-carrying qubit will be useful. In general, by treating one of the qubits as 
the data qubit (d) and the remaining ones as ancillae (ai, a2), the partial trace operation over the ancillae 

Pi7i = Trai.aslftn} ^ Tra^ (ftn ) } = Q{Ptn) (13) 

associates to a three-qubit process £ a one-qubit process Q on the data qubit d alone, provided that the latter is initially uncorre- 
cted with the ancillae and the initial state |ai a.2) is known. For the experimental realization, the choice |ai a.2) — |00) will be 
relevant. It is then readily seen that the above £z process corresponds, for instance, to applying a one-qubit map of the form 

ElpEl+EtpEt , (14) 

where E^ — ( 1 ± ct^ ) /2 are the usual z-idempotents, {E^ ) ^ = E^ . As expected from physical intuition, this process is nothing 
but crusher phase damping on the (un-encoded) data qubit. 



D. Non-abelian error models 



Whenever two non-commuting error generators are relevant, the interaction algebra describing the resulting error process 
is non-abelian. In practice, we shall be interested at error models obtainable by using abelian noise processes as building 
blocks. Let £x, £y, £z denote crusher dephasing about x, y, z, respectively, with associated error algebras Ax, Ay, Az, as 
above. A simple way for inducing non-abelian error processes is through the sequential composition of abelian errors along 
different axes. In particular, crusher isotropic collective decoherence corresponds to cascading crusher noise processes about 
all three axes [Q]. For instance, £y£z£x — £yzx, with Kraus operators specified by the 64 (linearly dependent) products 
K^K^K^, a,b,c = 0, ... ,3, of the operators given in the previous subsection. While composite noise processes of this sort 
may not naturally occur in physical systems, they can be enforced in liquid-state NMR using readily available non-unitary control 
methods to be described later (see also [p4|]). In terms of error algebras, the operators describing the composite noise process 
£yzx can be thought to belong to an error algebra Ayzx which arises from the multiplication of the single-axis algebras, 

Ayzx = AyAzAx = span{ if^^ if,^ , K^K^Ry} C Ac . (15) 

In fact, one can show that Ayzx is the full collective algebra Ac by checking that a basis Sa of permutation-invariant operators 
(such as the one given in (^)) is contained in Ayzx- It also turns out that the full Ac may be generated from the composition of 
abelian error processes involving two non-commuting axes, e.g. £z£x — £zx- 

By reasoning as in the abelian case, the one-qubit map Q describing the effect of a crusher composite noise processes at the 
single-qubit level can be derived by evaluating the partial trace (|l3|) on the appropriate sets of 16 (for two-axes noise) or 64 (for 
three-axes noise) Kraus operators. The result is effective full depolarization on the physical data qubit, corresponding to a map 
of the form 

u—x.y,z 
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E. The role of the self-Hamiltonian 



From a physical point of view, it is worth pointing out that while noise processes that involve both quantum decoherence 
(phase damping) and dissipation (amplitude damping) always correspond to non-abelian error algebras, a non-commutative 
error algebra does not necessarily indicate the presence of genuine energy dissipation in the system. In general, in the absence of 
an internal Hamiltonian, Hs = 0, every noise process is effectively equivalent to (adiabatic) decoherence. For instance, because 
all the basic operation elements in ( [l2| ) are Hermitian, crusher collective dephasing as defined above is, consistently, a unital 
process, describing damping of phase information in the z basis. Composition of dephasing operations along non-commuting 
axes results in depolarization, which is still unital and can be regarded as decoherence in all bases. 

Because Hs is rarely zero in real systems, the interplay between the internal dynamics and the actual error generators J a turns 
out to be crucial for characterizing the overall open-system dynamics. Clearly, it is always possible to construct the interaction 
algebra A by including Hs, as done so far, among the defining interaction operators. In particular, the assumption Hs — eJz 
is implicit in the original definition of the collective noise model [^. Regardless of whether Hs commutes or not with the Ja, 
if a non-trivial NS is supported by the algebra A constructed in this way, this NS is completely stable against time evolution. 
Therefore this procedure is, in principle, ideal for devising robust quantum memories. 

However, including Hs among the defining operators for A is not desirable when the A becomes irreducible (thus without 
useful symmetries) or when Hs is regarded as a resource for effecting noise-protected manipulations of information. Suppose 
that an error algebra A is defined starting from the 1 and the Ja alone. Then Hs may or may not belong to A'. From the 
point of view of reliable QIP, situations falling into the first category are the most favorable, as the natural Hamiltonian directly 
implements a non-trivial logical evolution on any NS supported by yl 0, If Hs ^ A' instead, then Hs may still 
preserve a given NS 0, but most likely it will have the undesired effect of causing leakage outside the intended space. 
While various schemes are available in principle to cope with these effects [|l], ^], it is not clear to what extent these may be 
viable with realistic control resources. 

In liquid-state NMR systems, collective error models may naturally play a role in describing relaxation from fully correlated 
fields in homo-nuclear species ^T^ . However, the assumption that Hs is proportional to is invalid due to the chemical 
shift effects. For the specific molecule used in the current implementation, Hs (to be explicitly given later) will turn out to 
satisfy none of the properties of belonging to the full collective error algebra Ac or its commutant A',.- Nevertheless, the system 
may be used to demonstrate how robustness against Ac can be achieved by constructing an appropriate NS. 



F. The three spin- 1/2 noiseless subsystem 

Because, as noted earlier, the eigenvalue j of the total angular momentum is conserved, simultaneous eigenstates of J^, 
are a natural basis to describe the state of the three particles. The possible values of j are j = 3/2, 1/2, corresponding to a 
decomposition of H as the direct sum of two invariant subspaces, H ~ 7^3/2 ffi 'H1/2, respectively. The quantum numbers 
j, iz suffice for completely labeling basis states in 7Y3/2, which is the four-dimensional subspace spanned by totally symmetric 
states {i.e., states in 7Y3/2 transform under particle permutations as the one-dimensional symmetric irreducible representation 
of S^). However, with j = 1/2 and jz = ±1/2 this is no longer true for the subspace Tii/2, which is also four-dimensional 
as the eigenvalue j = 1/2 is in fact doubly degenerate. Physically, this degeneracy accounts for the fact that there are two 
distinct paths for obtaining a total angular momentum j = 1/2 out of three elementary 1/2 spins. Let ^ = 0, 1 be an additional 
quantum number that labels these two possible paths. Because collective errors do not have access to the quantum numbers 
of the individual spins, and the resulting global quantum numbers are the same in both paths, the noise can neither distinguish 
which value of i is realized, nor change that value. Thus, £ corresponds to a conserved, two-dimensional DOF under the noise. 
In fact, this is the NS we are seeking [|[ ^ |4^ ]. 

More formally, basis states in 7ii/2 are labeled by two quantum numbers, \i, jz), with £ — 0,1, jz — ±1/2. For fixed £ (fixed 
path), the resulting subspace carries a copy of the two-dimensional irreducible representation of the angular momentum group 
su{2) corresponding to j = 1/2. There are two such copies, and operators in Ac do not mix them and act identically on both. 
For fixed jz, one obtains instead a copy of the two-dimensional irreducible representation [2 1] (in Young tableau notation) of the 
permutation group ^3 There are again two such copies, and now operators in A'c do not mix them and act identically 

on both [|50|]. With respect to the interaction algebra Ac, we can thus identify the DOF corresponding to £ with a subsystem L, 
which is fully protected against errors, and the DOF corresponding to jz with a syndrome spin- 1/2 subsystem Z, which instead 
experiences the errors. We write 7Yi/2 ^ Hl <^ H-z, and make this identification explicit through the following correspondence 
with states in the computational basis [Q|: 

|0)l(8)| + 1/2)z ^ -i=(|001)+cj|010)+cj2|i00)) , 

|0)L(g)|-l/2)z ^ 4^(|110)+cj|101)+cj2|011)) , 
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\1)l®\ + 1/2)z ^ -^(|001)+c^2|010)+Lj|100)) , 

|1)l®|-1/2>z ^ -^(|110)+c^2|101)+L^|011)) , (17) 

where lu — exp(«27r/3). Note that two of the totally symmetric states spanning 7^3/2 are |000), while the remaining 

two are obtained from ( jl7| ) by dropping the uj^uj'^ phase factors. Thus, the full set of basis states for 7^3/2 and 7ii/2 provides 
an explicit realization of the general state space decomposition given in (Q), with U being the rotation needed to bring the 
computational basis to the appropriate angular momentum basis, and = 8, mi = 1, di = 4, TO2 — 2, d2 — 2, respectively. 
Equivalently, the NS qubit living in TY^ can be identified by combining the 1 and the invariant operators ( [lO| ) into three operators 

ai^^ e A'^ that behave algebraically like the Pauli operators. This gives [|[] 

= \{^+ -12) A/2 , = ^(-^23 - S,,)P,„ , ) = ) , (18) 

where P1/2 — t/2— (si2 + S23 + S3i)/6 denotes the projector onto the subspace 7ii/2. Note that the logical observable 
is simply the restriction to 'H1/2 of the permutation 7ri2 swapping qubits 1 and 2. 



G. Collective operators in the noiseless subsystem/syndrome representation 



It is instructive to take a closer look at the action of the error generators J„ and of some relevant collective error processes 
directly in terms of the decomposition 

un^n3/2®ni/2^n3/2®nL<E)nz , (i9) 

where U effects the change of basis mentioned above and the subsystem identification within Ti.i/2 is explicitly given by (p^). 
By using the identity 1 + w + = 0, it is easy to verify that the restriction of the collective noise generators to the 7Yi/2 
subspace acts as follows: 

2J^Pi/2 ^ l(^) ® (-fxf , 

2J,Pi/2 <^ Il(^)®(+<Tf)). (20) 

Thus, the J„ act as single-qubit errors on the syndrome subsystem alone. Explicitly, this means that if, for instance, the initial 
state of the system is given by 

i'i/2l*>i23 = I*) - («|0) + ® (7I + 1/2) + - l/2))z , (21) 

for appropriate coefficients, then the result of a collective rotation by, say, dy about y is 

e-*^^«|^) ^ (a|0)+/3|l))i®e+*«»"«^'/2(^|^^/2)+<5|-l/2))z 

= (a|0) + (3\1))l ® ((7COs(0j,/2) + Jsin(0j,/2))| + 1/2) + (5 008(^^/2) ~~ jsm{9y/2))\ - 1/2))^ , (22) 

while a collective z rotation is simply 

^-^o..h\^^ = (a\0)+(3\l))L®e-'''^'"'^\j\ + l/2)+S\-l/2))z 

= (a|0) + /3|1))l ® (7e-^'^/2| + 1/2) + Se+''^/^\ - l/2))z , (23) 

and so on. Two observations are worth making. First, if either (5 or 7 is zero, then ( ^3| ) is a direct manifestation of the fact 
that each of the pairs of states in (|l7|) with fixed jz is a one-qubit DPS against pure z noise [|o|. As explained in Sect. IIIB, 
initialization of the ancillae qubits aslai 3.2) — |00) will correpond to encode into the = —1/2 subspace, thus 7 = in 
the implementation. Second, Eqs. ( p2[ ) and ( |23| ) together show that encoding into the Z subsystem would instead result in a 
qubit fully controllable in terms of homogeneous local unitaries [|T]| i.e., transformations of the form ?7'^^ (E) J/'^^ ® J/'^' on the 
physical qubits - for instance, non-selective ("hard") tt pulses about two non-commuting axes. 
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A procedure similar to the one outlined above may be applied to picture the effect of an arbitrary, known error model. Take, 
for instance, crusher collective z-dephasing with Kraus operators = {K^} given in (|l2|). One finds that 

ZZZPi/2 ^ l(^) (g) (-(t(^)) . (24) 
Thus, when restricted to 'H1/2, the action of both K§ and A'f is zero, whereas 

K2P1/2 ^ 1^^^<»E^_^^ , 

KIPi/2 ^ l'^^ . (25) 

This just means that the action of crusher collective phase errors on the physical system can be pictured as a crusher phase 
damping channel on the syndrome subsystem Z alone. Because, as we shall also comment later, the state of the latter abstract 
subsystem is mapped, upon decoding, into the state of a physical ancilla qubit carrying the error syndrome [|6[, this means 
that the physical syndrome subsystem will have experienced full phase damping under the same conditions. By the same 
reasoning, the action of crusher composite noise processes can be understood in terms of the composition of phase-damping 
channels affecting the Z subsystem along various axes, translating into depolarization of the physical syndrome subsystem. 
These observations will be corroborated by experiment. 



H. Verifying infinite-distance error-correcting properties 



As a result of the above analysis, quantum information encoded in the L subsystem is protected indefinitely in time, without 
requiring any active intervention. In the language of QEC [[TT|], this stability against the full collective algebra Ac characterizes L 
as an infinite-distance QEC code for arbitrary collective errors. Formally, this follows from the fact that to the NS one can asso- 
ciate a QEC in the usual (subspace) sense by choosing an initial reference state \e)z in the syndrome subsystem (corresponding 
to "no error"), and by letting the code subspace C be defined by 

C - span{|0)c = |0)l ® \e)z , |l>c = |1>l ® \e)z} ■ (26) 
Then the basis states of C verify the necessary and sufficient conditions for recovery from all errors in Ac [ |2^ |5^ , 

{ic\ElEb\jc) = aab6^J , VEa, EbeAc, (27) 

for appropriate coefficients aat - independent of the logical index Note that, a priori, infinite-distance behavior as expressed 
by ( p7| ) applies regardless of whether L is supported by a DPS or by a proper NS. For a DFS, the syndrome state is fixed thus it 
can be effectively disregarded. For a NS, however, the fact that errors are allowed to evolve the state \e)z non-trivially implies 
that the latter may be effectively arbitrary. In both situations, no recovery is needed for maintaining information in L |pl]]. 

In a NS-QEC experiment, one is interested at inferring robustness properties of the encoded information by looking at the 
code performance under a given set of quantum processes. Suppose that, as mentioned in the previous subsection, verification is 
constrained to having initialized the syndrome subsystem in the state | — l/2)z- What kind of conclusions can one draw? The 
analysis is relatively simple under the assumption of perfect implementation fidelity. The relevant points can be summarized as 
follows: 



Verifying that quantum information is preserved under the implementation of an error process with Kraus operators {Ka} 



implies stability under any error operator K e spanjii'a} 125 1 



• If the set of correctable error operators contains an error algebra A i.e., spanjA'a} 12 A, the implementation verifies 
infinite-distance behavior under A [pl]|. 

• Verifying a NS under A requires verifying infinite-distance QEC for every possible initial state of the syndrome subsystem. 

These observations can be applied to analyze both abelian and non-abelian collective error models. For instance, suppose 
we observe stability under crusher x noise. Then we can conclude that operators in Ax have identity action on L when Z is 
initialized to | — 1/2)2. However, because the state | + 1/2) ^ can be reached from | — 1/2) ^ by application of error operators 
in Ax, one has effectively verified a NS against Ax- While a similar argument applies to y noise, the state of the Z subsystem is 
preserved under z noise. Thus, observation of stability under Az for fixed initialization in \ — l/2)z only implies the verification 
of a DFS-behavior under Az. A proper NS-behavior under the full collective Ac can be inferred, in principle, in various ways. 
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Keeping the experimental | — 1/2) z preparation constraint in mind, the simplest procedure is to ensure stability under a family 
of quantum processes whose sets of Kraus operators globally span Ac- For instance, one can check that 

span{K:K^ , K^K^ | a, 5 = 0, . . . , 3} = ^ (28) 

for any choice of a composite crusher process which involves two non-commuting axes ii, ii e.g., E^y and Eyz- The implemented 
set of processes are discussed in Sect. IVB-C. One may notice that the error operators describing Eyz are obtainable as a subset 
of the errors operators induced by £yzx- Thus, the implemented set is sufficient to infer that a NS for the most general collective 
noise has verified, at least in the limit of sufficiently high fidelity. Establishing what minimum fidelity threshold is required for 
inferring NS-verification under non-ideal conditions is a separate interesting issue, whose analysis is beyond the scope of this 
work. 

The three-qubit NS turns out to provide the smallest code capable to correct one qubit against the full Ac- The same, infinite- 
distance protection can be accomplished by using a DFS, but the most efficient DPS requires four physical qubits [Ql. It is 
worth noting that, for a given physical system, more efficient codes may exist if additional symmetries are present beside the 
permutational one. If, for instance, axial symmetry also applies (i.e., the error model belongs to the class of collective single-axis 
phase damping), then a three-dimensional subspace of Ti, may be protected with infinite distance by a DFS (corresponding to a 
fixed jz eigenvalue, jz = ±1/2). The situation is summarized in Table I. 



III. EXPERIMENTAL DESIGN AND METHODS 

A. Liquid-state nuclear magnetic resonance 

All experiments were performed using liquid-state NMR techniques on a sample of ^■^C labeled alanine (Fig. |l]) in D2O 
solution, using a 300 MHz Bruker AVANCE spectrometer. The spin's evolution is governed by the internal Hamiltonian Hs, 
which in the weak coupling limit is accurately described by 

3 

fe=l,2,3 j>k=l 

where Vk represents the chemical shift frequency of the fcth spin, and Jkj the coupling constant between spins k and j, respec- 
tively (in Hz units). Radio frequency (RF) pulses are used to modulate the dynamics to produce the desired net evolution in 
the spin frame where the above internal Hamiltonian is defined. The interaction with the control field generated from a single 
transmitter has the form 

Hc,t{LORF, 0, = E e-^("«-*+*)^^'/2(^ai72)e^("«-*+«'^^'/2 , (30) 

k 

where the transmitter's angular frequency LUpp, the initial phase 0, and the power uo are tunable over an appropriate parameter 
range. NMR QIP has been extensively discussed in the literature [^. At room temperature, NMR qubits exist in highly 
mixed, separable states and so NMR QIP relies on "pseudo-pure" (p.p.) states whose traceless (or "deviation") component is 
proportional to that of the corresponding pure state. The identity component of the density matrix is unobservable and, if the 
dynamics are unital, constant. The assumption of unital behavior has been validated experimentally, see Sect. IVD. Under these 
circumstances, the evolution of a p.p. state is equivalent to the one of the corresponding pure state. The 3-spin p.p. ground state 
Qz = 1 000) (000 1 was generated using standard gradient-pulse techniques whose details can be found in [|4p. State preparation 
was experimentally verified by tomographically reconstructing the 3-spin density matrix [|5^. A constant amount of identity 
component was added to all reconstructed density matrices such that the ground state fidelity with respect to the intended 3-spin 
p.p. state was optimized. 



B. Encoding and decoding quantum networks 

It is essential that the experimental procedure is designed to allow for the protection of an arbitrary one-qubit state, \4'm) = 
ajO) + /3|1), with a and /3 potentially unknown. Because collective errors affect the Z subsystem hence induce a non-trivial 
evolution the encoded states, the decoding transformation Udec must map the entire set of encoded basis states ( [T^ back to 
the computational basis properly. Thus, a good decoding transformation provides an explicit realization of the mapping U^^ 



appearing in (|19[). Given Udec, an encoding transformation may be obtained by letting Uenc — U^^\. Various choices are possible 
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in principle, differing in the identification they establish between the state space of the abstract L and Z qubits and the physical 
qubits. The choice of Udec we implemented maps L, Z into the data qubit 2 and the ancilla qubit 3, respectively: 



0)l ® 


| + l/2)z^ 


|J 


= l/2,^ = 0,j, = 


+ 1/2) H 


-> |001) 




1)l ® 


| + l/2)z^ 


|J 


= 1/2,^= - 


+ 1/2) H 


- 1011) 




0)l ® 


|-l/2>z- 


|J 


= l/2,^ = 0,j, = 


-1/2) H 


-> |000) 






|-l/2>z- 


|J 


= 1/2,^= - 


-1/2) H 


-> |010) . 


(31) 



The fact that (^T]) only specifies Udec on the 7i 1/2 subspace is reflected by the fact that the value of the first ancilla qubit remains 
set to zero. Udec is uniquely determined by also defining its action on the \ j — 3/2, j^) states spanning 7^3/2 or, equivalently, 
the mapping with the remaining four computational basis states: 



b--3/2,j, = 


+3/2) h- 


|100) 




b- = 3/2,j, = 


+ 1/2) ^ 


^ -i|lll) 




b = 3/2,j, - 


-1/2) ^ 


^ -i|110) 




b = 3/2,j, = 


-3/2) ^ 


. |101). 


(32) 



Having determined Udec, no general efficient procedure is known for designing a logical network of realizable one- and two- 
qubit operations effecting Udec- The implemented networks for Udec and Uenc are shown in Fig. |[ Due to the weak strength of 
the Ji3 coupling, no direct gates between spins 1 and 3 are used. For practical convenience, Uenc is simplified by using qubit 3 
as the information carrier (i.e., switching qubits 2 and 3 with respect to the output) and by taking advantage of the knowledge of 
the starting state |00) of the input ancillae. Thus, the implemented Uenc effects the transformation 

Uenc\m^n) ^ («|0)l ® \ - l/2)z + /3|1)l ® | - l/2)z) = \i^^n) L ® \ ' l/2)z • (33) 

The logical gates involved in the NS encoding and decoding circuits were first mapped into ideal pulse sequences via standard 
quantum network methods [|3]|. Pulses were then implemented by strongly modulating the internal Hamiltonian (^9|) of the 
alanine molecule with externally controlled RF magnetic fields as mentioned above. Details of pulse design can be found 
in [^4]]. Combinations of rotations that were used in multiple places in the sequence were merged into a single pulse and 
directly implemented. In practice, inaccurate preparation of the ancilla state |aia2) = |00) and operational errors may result 
in unintentionally populating states in the 7^3/2 subspace. Given the decoding action (P2|), this could contribute to observable 
signal upon discarding qubits 1 and 3. The fact that contributions originating from 7^3/2 remained negligibly small was inferred 
in the implementation from the absence of appreciable double- and triple-coherence decay modes and from the stability of the 
observed signal against applied noise strength [ pO| (See also Sect. IVB). 

It is worth emphasizing the difference between the ability to encode an arbitrary quantum state and the ability to accomplish a 
desired state preparation. While the latter is appropriate in the context of initializing a quantum algorithm, the former is crucial 
for quantum memory purposes (see also p9|]) and for QEC in general. In the specific case, much simpler procedures would 
suffice for initializing the system into a desired, known NS state. For instance, one can observe that with an appropriate choice 
of basis a NS always contains a state which is the tensor product between a two-qubit singlet and a one-qubit logical state (see [|[] 
for explicit encoded states alternative to ([l7|)). Then, similar to the case of initialization in a DFS [|56|], synthesizing a NS state 
may be achieved by relying on a unitary transformation that prepares singlet states. Alternatively, one could exploit appropriate 
non-unitary control such as cooling, as suggested in [^. For either DFSs or NSs, the existence of such initialization procedures 
does not automatically translate into the existence of an efficient way for effecting a general encoding. From this point of view, 
a systematic comparison of network complexity for DFSs vs NSs is worth being examined in more detail, and will be addressed 
elsewhere. 



C. Heisenberg representation of collective noise 

As discussed in Sect. IIA, it is unlikely that noise symmetries directly imply the preservation of one (or more) of the natural 
subsystem's DOF (such as the states of the physical spins). Yet, an abstract information-carrying subsystem not identifiable 
with any of the physical qubits exists through the subsystem identification (|l7|) 26 1. The decoding operation corresponds 



to extracting the abstract protected DOF by mapping it to the natural DOF associated with the data qubit. In the abstract 
NS/syndrome representation, the desired error-correcting behavior translates into the property that collective errors act only on 
the Z subsystem, as discussed earlier and explicitly verified in (|o|). An equivalent description can be constructed directly in 
terms of the physical DOFs by looking at the error algebra Ac in an appropriate Heisenberg representation determined by Udec- 
Let Ea e Ac be any collective error operator. The state of the entire system after it has been encoded into the NS, affected by 
the collective error, and then decoded from the NS is given by 

I V'o«t )d 1 0)ai I (/'out )a2 = t^dec-Ea[/]eclV'm)d|0)ai|(/'m)a2 = \'^in) d\Q) B.A4>in) ^-2 , (34) 
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where the Heisenberg-transformed error operator is defined by = UdecEaU\^c- Because Ea is in Ac, the action of E^ 
can be inferred by knowing the transformed collective generators J^, u — x,y, z. Using the decoding network of Fig. ^, and 
recalling that d, ai, a.2 denote the data, the first and second ancilla qubit, respectively, one finds 

2jf = i?;n-'^?^)l' + S-^'^?(l'-2cos(V3)^^-2sin(V3)^,') 

2jf = S;i(-CT^2)ll^ + S^V^^(]l^-2cos(V3)cT^-2sin(V3)cTf) 

2jf = El'(+al^)i'' + E^Jal^ii" + 2ai) . (35) 

The above equations make it explicit that the action of any collective error is identity on the data bit, provided that the state of 
the first ancilla ai qubit is properly set to |0). This identity action is the counterpart of the one expressed in the Schrodinger 
picture by Eqs. (|20|). On the other hand, the second ancilla qubit a.2 is in general evolved non-trivially by the noise. According 
to (^5]), the protection afforded by the NS can be equivatently understood as a rotation of the error algebra such that the abstract 
protected DOF is identified with the data qubit conditionally to the state of the first ancilla qubit. 



D. Design of the evolution period 

In order to verify the behavior of the NS code in a controlled setting, the delay period between encoding and decoding needs 
to be carefully designed so as to implement an effective evolution under a precisely known error model. An important building 
block for the procedure is the ability to implement a "no-op" evolution i.e., an evolution corresponding to the identity operation 
that provides the reference of no applied errors. 

Additional constraints exist because, as mentioned, the assumptions made on Hs in Sec. ^ are not achievable in the current 
experimental setting. In particular, Hs does not itself respect the permutational symmetry of the collective error model, neither 
does it belong to the commutant A'^ of the collective error algebra, causing the system to depart from the protected NS. The 
net, unwanted evolution induced by Hs over the delay period can be averaged out by using refocusing techniques [ p4| , The 
implemented sequences of tt pulses are depicted in Fig. ||. 

Because the system evolves through intermediate states outside the NS during the refocusing cycle, protection of the quantum 
data cannot be expected if noise is acting over the entire evolution period. However, what is crucial for verifying the robustness 
of the intended NS is that the information resides inside the protected space while the noise is applied. This can be achieved by 
making sure that the basic building block for engineering collective noise processes is a z-gradient, in which case Hs € A'^, and 
by applying the noise only during the sub-interval of the decoupling cycle corresponding to the identity frame (i.e., to the free 
evolution sub-interval) [|l7, 41,5^. While this suffices for testing collective dephasing processes along z, error algebrs of the 



form Ax, Ay can also be probed, by simply sandwiching a z-noise process with the appropriate collective rotation pulse about 
either the x or y axis Finally, cascades of noise blocks involving different axes were implemented by applying multiple 
noise blocks in series during a longer evolution period. 



E. Implementation of collective error models 

Two different engineered noise models were implemented to test both the weak and strong noise limits [|^. Both relied on 
linear magnetic field gradients in order to create an incoherent evolution over the spatial distribution of the sample [|60|]. The net 
phase evolution caused by a gradient pulse for a spin located 5z from the center of the sample is given by 

A$(,5,) = exp ( - «7^'5,) , (36) 

where 7 is the gyro-magnetic ratio of the nuclear species and determines how strongly the magnetic moment of the spin couples 
to the magnetic field. A strong gradient pulse causes spins at the edge of the sample to evolve through many cycles, producing an 
almost uniform distribution of phases across the sample. In the absence of molecular motion, the effects of this incoherent evo- 
lution could be refocused by a reverse gradient pulse. However, in liquid samples at room temperature, the spins are undergoing 
a random spatial diffusion that renders this incoherent process effectively irreversible, resulting in unrecoverable loss {i.e. deco- 
herence) of the quantum information. The exponential attenuation of the signal associated with the combined gradient-diffusion 
process is described by a factor [|o[, 

A = exp - £> j k^{t)dt^ , (37) 

where Z? is the molecular diffusion coefficient and = ^t{dBz/dz). For a gradient pulse of duration 5, and a diffusion period 
of duration A followed by an inverse gradient also of duration 5, the total attenuation of the coherence reduces to 

A = exp ( - Dj^{dB,/dzfS^{A + —)j . (38) 
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Therefore, one can associate this evolution with a r2-like process, exp(— r), where 1 /r represents the relevant, effective noise 
strength. By equating A = exp(— t/r), the elapsed time t — A + 26 being the duration of the gradient-diffusion process, the 
noise strength is given by 



2<5 



i = Dj^idBjdz fS' . (39) 

This reduces to the simpler expression quoted in [^l]] if 25 <C A. According to (^9|), the noise strength can be tuned by 
changing either the holding duration A or the gradient strength dBz/dz. In order to keep the effects of other natural decoherence 
mechanisms fixed in different experimental runs, only the gradient strength was varied in practice (decoherent implementation). 
Each effective noise rate l/r was independently measured using a stimulated echo sequence [|60|]. 

If no reverse gradient is applied to refocus the magnetization, then the evolution remains incoherent until the ensemble signal is 
acquired, at which time the spatial degrees of freedom are traced over rendering the evolution effectively irreversible. Therefore, 
the crusher noise limit can be probed by using a strong gradient pulse with no reverse gradient pulse (incoherent implementation). 
The details of the no-op sequences and the different noise implementations are given in Fig. ^ 



IV. RESULTS 



A. Metric of control 

We used the entanglement fidelity [|2|| as the reliability measure quantifying how well quantum information was preserved 
under the evolution of an implemented super-operator Q. Given an operator-sum decomposition of Q in terms of Kraus operators 
{Ap}, and an input state p, the entanglement fidelity can be computed as 

Fe{Q,p) =Y.^r{pA^)Tr{pAl) . (40) 

For a uniform distribution of input states, the maximally mixed input state, p — 1 /2, is used to characterize the un-biased channel 
performance. The entanglement fidelity can then be inferred from experimentally available data once it is expressed in terms of 
either pure state fidelities or relative input-output spin polarizations. Under the additional assumption of unital dynamics, i.e., 
= 1, the relevant expressions for single-qubit quantum process tomography are, respectively: 

Fe{Q.) = \(F\+,) + F\+y) + - l) , (41) 

or, equivalently, 

Fe{Q) = ^[l+Px+Py+P.) ■ (42) 

In (plj), the pure state input-output fidelity Fy^^ = Tr{\ijj) (iplQdip) (ipl)}, and \u), u ~ x,y,z, are the eigenvectors of the 
corresponding Pauli operator with positive eigenvalue (thus, 1 + 2) =J0), \ — z) = |1)) ^sj]. In (p^, p^ — Tr{(T„Q((T„)}/2 



represents the relative output polarization given input au [M, 61 63]. The implementation allowed to explicitly validate the 



consistency of the above expressions within experimental accuracy. 



B. Weak collective noise along a single axis 

The ability of the NS to protect quantum information against weak collective noise was tested using the gradient-diffusion 
techniques described above. The entanglement fidelity of the data qubit was experimentally determined from Eq. ( |4l| ) as a 
function of noise strength for single axis collective noise. Separate experiments were carried over for extracting the behavior of 
both the NS-encoded qubit and the un-encoded data spin. The results are summarized in Fig. Q The measured curves are fitted 
to a decaying exponential function, Fe = A exp(— t/r) + B, as expected for a dephasing channel induced by a single-axis noise 
mechanism. The un-encoded data's fit is characterized hy A — 0.51 ± 0.04 and B — 0.43 ± 0.03, confirming the expected 
decaying contribution. The NS data's fit is instead characterized hy A = 0.03 ± 0.03 for both y and z noise, with constant 
coefficients B = 0.62 ± 0.02 and B = 0.64 ± 0.02, respectively. 

For all implementations, departures from the expected ideal behavior may be explained by pulse imperfections as well as by 
naturally occurring relaxation processes, whose effects are assumed to be independent of the applied noise strength. First, the 
fact that the decaying contribution (coefficient A) is small for both of the encoded situations indicates that the data resides in the 
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NS during the application of the noise, which makes it insensitive to collective errors. Second, the fact that the constant term 
(coefficient B) is greater than 0.50 [Q, confirms that quantum information is retained, in principle, for arbitrary noise strengths 
despite significant imperfections in the implementation of the encoding and decoding operations. Finally, the extrapolation of 
the data bounds the performance of the NS in the strong noise limit. 

As discussed in Sect. IIG-H, unlike the case of a noiseless subspace the encoded NS states are not eigenstates of all noise 
operators. This implies, in general, a decay of the full-state fidelity due to evolution of the ancillae qubits under the action of 
the noise. Because the system always resides, within the experimental accuracy, in the spin- 1/2 subspace, the ancilla ai remains 
unchanged by the application of the noise while the ancilla qubit &2, which is mapped from the syndrome Z, is decohered by the 
action of y noise operators. For the case of collective z noise, and with initialization in a fixed jz subspace, the system resides 
into a DFS, as noticed repeatedly. Thus, the state of the second ancilla qubit remains also unchanged under the z-dephasing 
process. These behaviors are experimentally confirmed in Fig. || Note that because |a2) is always initialized to |0), does not 
provide an appropriate metric. An average input-output fidelity is evaluated instead, resulting from a uniform average over the 
data input states. 

C. Incoherent implementation to mimic strong collective noise 

A variety of incoherent collective noise processes were also implemented to explore the strong noise limit and to establish 
robustness of the implemented NS against the full, non-abelian Ac- The entanglement fidehty data, again calculated via Eq. 

are presented in Table ^ [^0|]. As expected, the un-encoded qubit's entanglement fidelity drops to 0.50 for single axis noise 
(corresponding to a full-strength dephasing channel), and to 0.25 for a cascade of non-commuting noise blocks (inducing a full- 
strength depolarizing channel). For the NS-encoded case, the entanglement fidelity again departs from unity due to pulse errors 
and non-collective natural decoherence mechanisms, but it is changed only slightly under the action of different noise blocks. 
The significant improvement in the amount of information retained under a cascade of strong noise mechanisms confirms the 
expected benefits of the NS encoding already indicated by the weak noise data. 

D. Experimental determination of channel super-operators 

To gain further insight about the quantum processes realized in the experiment, the direct reconstruction of the super-operators 
describing some representative crusher dephasing channels was also obtained, leading to explicit sets of single-qubit Kraus 
operators. The measurement of the input-output relations for the set of input states | + x) {+x\, \ + y) (+y|, | + z) {+z\, and | — 
z) (— z| is sufficient to allow the determination of a Kraus form for a desired super-operator [pSl pq]. While this analysis assumes 
that the channel is trace-preserving, unitality is now not assumed, thereby allowing to independently test this assumption. 

Table [ll| collects the experimentally determined one-qubit super-operators and corresponding entanglement fidelities, calcu- 
lated via Eq. ( p^ (p = 1/2) for five representative channels. Because the diagonal elements in each of the channel super- 
operators measure the relative polarizations p„ introduced in (^2|), Fe{Q) can equivalently be evaluated from the trace of the 
relevant super-operator representation (up to a factor 1/4). 

First, the entanglement fidelity values are consistent with the ones calculated via Eq. (p4|). Second, the data show that the 
unitality assumption is broken by an average deviation no larger than 5%. Finally, it is worth stressing that the measured super- 
operators contains complete information on the relative contributions of coherent vs decoherent errors. For instance, the purity 
of an output state pout, given by ^ = T^^iPout} ^ given input state 1 , , cr^ , cr^ can be inferred from the sum of the squares 
of the elements of each row of the super-operator 

V. CONCLUSIONS 

We presented a thorough theoretical and experimental investigation of the significance of a NS in the context of quantum 
information protection. This exploration demonstrates both the utility of the NS and its limitations. In particular, the assumption 
that the internal Hamiltonian is either proportional to one of the noise operators or that it is trivially zero is not valid for the NMR 
implementation, and is unlikely to be met exactly in any experimental realization of a quantum information processor If such 
assumptions are at least approximately met, NSs represent the most efficient means for protecting quantum information in the 
presence of noise interactions with dominant symmetry components. For the example of the 3-qubit collective noise model, our 
experimental results convincingly demonstrate improvements in protecting information against a class of both abelian and non- 
abelian collective error models. Finally, experimentally reconstructed single-qubit super-operators provide detailed information 
about the properties and the relative contributions of coherent and decoherent errors in the implementation. 
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Axial independent 
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Axial independent 
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3bit QEC 
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General collective 


Arbitrary 
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3bit NS 
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General collective 
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3bit NS 
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Axial collective 


Arbitrary 


Perm, axial 


4 


3bit DPS 


3 


Axial collective 


Weak 


Perm, axial 


2 


3bit DPS 


3 



TABLE I: Comparison between relevant error models on 3 qubits and corresponding error control strategies. The various columns list, in the 
order: the type of environment interaction (independent or collective); the relevant noise strength (or, equivalently, the order in time at which 
protection is sought (jll|], weak strength corresponding to first-order effects); the type of symmetry (permutation or axial); the dimension of 
an error basis for the appropriate error set; the most efficient error control code available for implementation on 3 qubits; the dimension of the 
state space C corresponding to the protected DOE. 3bit QEC means the standard QEC code for weak phase noise m. Note that finite-distance 
QEC could also be used to protect one qubit against general weak noise (either independent or collective), but the smallest such code requires 
5 qubits 




PIG. 1 : Molecular structure of C labeled alanine. The three carbon spins are used as qubits. In a reference frame that rotates at a frequency 
75.4736434 MHz, the internal Hamiltonian is accurately described by Hs = 7r[i/i(ji^' + :^20-i^' + z^acri^^ + {Ji2a-Pai^^ + Jizaf' af^ + 
Ji3cri^Vf ^)/2)], where i^i - i^o = 7167 Hz, V2 ~ = -2286.5 Hz, vz ~ ^ -4881.4 Hz, J12 = 54.1 Hz, J23 = 35.0 Hz, and 
J13 = —1.3 Hz. As it is impractical to directly utilize the coupling between qubits 1 and 3 due to its weak strength, conditional gates involving 
the pair of qubits (1, 3) were effectively replaced by sequences of operations involving pairs (1, 2) and (2, 3) in designing and executing 
quantum networks. The T\ relaxation times for the three spins are approximately 21, 2.5, and 1.6 s, while the T2 times are about 550, 420, and 
800 ms, respectively. 
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FIG. 2: Logical quantum network for the NS experiment. The information is initially stored in qubit 3, while the ancillae qubits 1 and 2 are 
initialized in the state |0). First, an encoding sequence {Uenc) is applied in order to map the initial state of the system into the NS, according 
to (^^. Next, the qubits are stored in memory while different noise processes are applied. Finally, the information is transferred to qubit 2 by 
the decoding sequence (Udec)- As noted in the text, Uenc is a simplified version of [/^^ obtained by exploiting the knowledge of the state of 
the ancillae qubits. 
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Qz.ns 


0.86 


0.72 


0.76 


0.67 


Qzx, un 
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Q000,ns 
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TABLE II: Experimental data for implementation of various collective error models in the strong noise limit. The first column lists the one-bit 
quantum channels realized in the experiment. In addition to the applied error model, £x,£y,Sz, £zx, £zy, £yzx, the channel label specifies 
whether (ns) or not (un) encoding and decoding procedures were implemented. The processes Qo.ns, Qoo.ns, Qooo.ns differ in the duration 
over which the net identity evolution is applied. For each process, the input-output fidelities f involved in the process tomography as well 
as the resulting entanglement fidelities Fe are listed. Statistical uncertainties are ~ 2%, arising from errors in the tomographic density matrix 
reconstruction. 
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FIG. 3: No-op pulse sequences for both weak and strong noise applied to the data qubit (d) and the two ancillae qubits (ai and 32). Black 
(white) boxes represent tt pulses about the x (— £) axis. The experimental delay times were 6 ~ 5.5 ms and t ~ 0.5 ms. The bottom line 
(g) denotes the magnetic field gradient sequence used to introduce either weak (top) or strong (bottom) noise. As explained in Sect. IIID, all 
gradients were along the z axis and were only applied during the free evolution sub-interval of the whole evolution period. Collective noise 
along other axis was implemented by collective rotations sandwiching the entire noise block. Cascades of noise along non-conmiuting axes 
were implemented by applying multiple noise blocks in series. 
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FIG. 4: Measured entanglement fidelities for single-axis collective noise in the weak noise regime along either the y axis (NS-encoded 
(squares) and un-encoded data (triangles)), or the z axis (NS-encoded data only (circles)). The decay of the un-encoded spin, Cs, was obtained 
by turning off the encoding and decoding sequences and by subjecting it to the noise sequence alone. The data are fit to an exponential decay 
function, with the interpolated (solid) and extrapolated (dashed) lines shown in the plot. 
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FIG. 5: Measured average state fidelity for the second ancilla qubit, a2, for single-axis collective noise in the weak regime. Because the 
encoding transformation Uenc maps the initial state into an eigenstate of (see (133)), the ancilla qubit's state is unchanged by the application 
of noise processes along z. This demonstrates that the encoding also maps the system into a one-qubit DFS against collective jz-dephasing 
|p(J]. The case of collective y noise demonstrates, however, that information is not stored in an eigenstate of a generic collective noise operator. 



21 



Q:c,un Fe = 0.48 

1 ax Oy (Tz 

/ 1.00 0.01 0.00 -0.01 \ 1 

0.00 0.92 0.06 -0.19 

0.00 -0.02 -0.02 0.02 Oy 

\ 0.00 0.07 0.00 0.02 / 



Q0,ns = 0.69 

1 CJx Uy Oz 

I 1.00 -0.03 0.03 0.02 \ 1 

0.00 0.53 0.29 -0.16 

0.00 -0.25 0.56 0.05 Oy 

\ 0.00 0.10 -0.07 0.67 / cr^ 



Qx,ns = 0.67 

1 (Jx CTy CFz 

I 1.00 -0.03 0.02 0.02 \ 1 

0.00 0.54 0.28 -0.21 cr^ 

0.00 -0.20 0.56 0.09 Oy 

\ 0.00 0.21 -0.07 0.58 / 



Qy,ns Fe = 0.69 

1 Ox Oy Oz 

/ 1.00 0.02 0.05 0.03 \ 1 

0.00 0.55 0.27 -0.16 Ox 

0.00 -0.18 0.61 0.11 Oy 

\ 0.00 0.16 -0.09 0.60 / cr^ 



Qz,ns Fe = 0.70 

1 Ox Oy Oz 

I 1.00 -0.10 0.00 0.03 \ 1 

0.00 0.56 0.29 -0.10 Ox 

0.00 -0.19 0.54 0.03 Oy 

\ 0.00 0.12 -0.07 0.67 / 



TABLE III: Experimentally reconstructed super-operators for the one-qubit channels. A fixed row represents the decomposition along the set 
of Pauli operators {t,Ox,Oy,Oz} of the output corresponding to a given operator t,Ox,Oy,o'[ o^ taken as input. Because trace-preservation 
is assumed, the first entry in each row is constrained to unity, and the remaining entries in the first column to zero. Non-zero off-diagonal 
elements on the first row indicate the non-unitality of the channel. Entanglement fidelities, calculated via Eq. ( p^ are listed on top for each 
channel, and are also consistent with the values calculated from Eq. ( |42[) . which equals Tr{Q}/4. The above data may be interpreted as 
providing a matrix representation of the underlying super-operator as in J66|] upon appropriately transposing rows and colums. 



